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ABSTRACT 

The contribution of stars in galaxies to cosmic reionisation depends on the star formation 
history in the Universe, the abundance of galaxies during reionisation, the escape fraction of 
ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the 
star formation rate and clumping factor during reionisation in a cosmological volume using a 
high-resolution hydrodynamical simulation. We post-process the output with detailed radia- 
tive transfer simulations to compute the escape fraction of ionising photons. Together, this 
gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. 
The strong mass and redshift dependence of the escape fraction indicates that reionisation 
occurred between z = 15 and z = 10 and was mainly driven by proto-galaxies forming in 
dark-matter haloes with masses between 10 7 M Q and 10 8 M©. More massive galaxies that 
are rare at these redshifts and have significantly lower escape fractions contribute less pho- 
tons to the reionisation process than the more-abundant low-mass galaxies. Star formation in 
the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these 
proto-galaxies only contribute when the part of the Universe they live in is still neutral. After 
z ~ 10, massive galaxies become more abundant and provide most of the ionising pho- 
tons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent 
enough to have a major contribution to reionisation. Although the stellar component of the 
proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be 
detected by the James Webb Space Telescope (JWST), these sources are brightest in the Ha 
and Ly-a recombination lines, which will likely be detected by JWST in deep surveys. 
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1 INTRODUCTION 

One of the major challenges in modern cosmology is identifying 
the nature of the sources responsible for reionising the Universe. 
The Gunn-Peterson trough in the spectra of high-redshift quasars 
indicates that the IGM was highly ionised at z < 6 ( |Fan et aT] 
2006). The integrated Thomson (electron) scattering optical depth 
to the surface of last scattering, r c , suggests that reionisation was 
well underway by z ~ 10.5 (Komats u et al.|2011| >. It remains un- 
certain which sources transformed the IGM into its highly ionised 
state. The most likely candidates are stars in galaxies ( [Faucher-] 
|Giguere et al.|2008| >. 

The contribution of galaxies to reionisation depends critically 
on the star formation history in the Universe, the abundance of 
reionising sources, the clumping factor of the IGM and the fraction 
of ionising photons that escapes into the IGM, the so-called escape 
fraction. High-redshift surveys that are probing the star formation 
rate of galaxies at 6 < z < 8 show that an escape fraction of more 
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than ~ 30% is needed for the observed galaxy population to pro- 
duce enough photons to keep the Universe ionised ( Robertson et al. 
|2010|[Bouwens et al.|2 012). Although there is some evidence for 
a redshift evolution of the escape fraction, such high escape frac- 
tions are not observed in the local Universe (Bland-Hawthorn & 
|Maloney|1999||Deharveng et al.|200T]|Heckman et al.|200T) and at 
higher redshifts jlwata et al ;2009 , Nestor et al.|201 l||Vanzella et al.| 
|2012| >. All the massive galaxies targeted by observations have av- 
erage escape fractions less than 20%, although samples are small. 
With this escape fraction, the ionising emissivity of the observed 
galaxy population is insufficient to maintain reionisation. However, 
the bulk of the star formation rate at these high redshifts likely 
arises in galaxies below the detection limit of current observational 
facilities jTanvir et al.||2012) . If the escape fraction evolves with 
redshift, these galaxies could provide the majority of photons for 
reionisation (Trenti et al.|2010||Kuhlen & Faucher-Giguere|2012J . 

Numerical studies ( |Gnedin et all|2008| |Wise & Ceiil|2009l 
|Razoumov~ & Sommer-Larsen 2010; Yajima et al. 2010) find es- 
cape fractions between ~ — 1, with possible redshift or mass 
dependence. In part the large differences between studies may be 



© 2012 RAS 



2 J. -P. Paardekooper et al. 



caused by numerical issues, because the radiative transfer simu- 
lations are computationally challenging, making approximations 
necessary. However, most studies targeted only a few objects and 
different studies focussed on different mass galaxies at different 
redshifts, making comparison difficult. In simulations of idealised, 
isolated galaxies Paardekoo per et aL| ( |2011} found the main con- 
straint on the escape fraction to be the dense gas in the star-forming 
regions, which provides an explanation for the large spread in es- 
cape fractions reported in previous studies. In addition, they found 
that the escape fraction can vary over several orders of magnitude 
over the lifetime of the galaxy, making it necessary to determine 
the contribution of galaxies to reionisation not only sampling the 
mass function, but also over a wide range of redshifts. 

In this letter we present results on the escape fraction and ion- 
ising emissivity from a large statistical sample of proto-galaxies 
in a high-resolution cosmological, hydrodynamical simulation. In 
combination with the derived clumping factor of the IGM in the 
simulated volume we assess self-consistently the contribution of 
stars in galaxies to cosmic reionisation. 



2 METHOD 

In the standard cold dark matter paradigm, at the relevant red- 
shifts for reionisation most ionising radiation is produced by stars 
in proto-galaxies forming in dark matter haloes of < 10 9 Mq 
(Bark ana & Loeb|200l||Choudhury et al.|2008) . We compute the 
ionising photon production and escape fraction in proto-galaxies 
in this mass range using the First Billion Years (FiBY) simulation 
suite (Khochfar et al. in prep.; Dalla Vecchia et al. in prep.). The 
simulation we use contains 2 x 684 3 dark matter and gas particles 
in a comoving volume of 4 Mpc on the side, with a gas-particle 
mass of 1250 M©. At redshift 6, the simulation reproduces the ob- 
served mass function of galaxies and star formation rates (Khochfar 
et al. in prep.). 

For the FiBY simulation we use a modified version of the 



OWLS code iSchaye et al.J 2010 >. Star formation follows a pres- 



sure law tSchaye & Dalla Vecchia 20081, where we assume pop- 
ulation (Pop) III stars form at metallicities Z < 10~ 4 Z (with 
Zq = 0.02) and Pop II stars at higher metallicities. Supernova 
feedback is modelled by injecting thermal energy that is efficiently 
converted into kinetic energy without the need to turn off radia- 
tive cooling temporarily (Dalla Vecchia & Schaye 2012 1. Feedback 
from reionisation is modelled as a uniform UV-background fol- 
lowing Haardt & Madau (2001) by switching from collisional to 
photo-ionisation equilibrium cooling tables. Gas above a density 
threshold of n s hi c id = 0.01 cm -3 is modelled to shield against 
ionising radiation (Nagamine et al. 2011). We assume that reioni- 
sation takes place within the bounds set by WMAP (Komatsu et al. 
|201 starting around redshift 12 and ending around redshift 9. 
This is consistent with the computations of the ionising emissiv- 
ity from proto-galaxies that we present in this letter. We will show 
that between redshift 12 and 9 the proto-galaxies in the simulation 
produce enough ionising photons to reionise the computational vol- 
ume. 

We have extracted all haloes from this simulation that contain 
at least 1 star, 1000 dark matter particles and 100 gas particles for 
post-processing with radiative transfer simulations. This results in 
more than 11000 haloes between z — 20 and z — 6. We deter- 
mine the escape fraction with an updated version of the SIMPLEX 
radiative transfer code (Paardekooper et al. 2010), computing the 
absorption of the ionising radiation by both hydrogen and helium 



atoms in 10 frequency bins until the photons reach the virial ra- 
dius of the halo. We will discuss these simulations in more detail in 
a forthcoming paper (Paardekooper et al. in prep). The luminosity 
and spectra of the star particles are computed from stellar synthesis 
models for both Pop III ( |Raiter et al.||2010"| > and Pop II ( |Bruzual| 
|& Cha riot 2003 1 stars. Combined with the escape fraction, we ob- 
tain the number of ionising photons that every proto-galaxy in our 
simulation contributes to cosmic reionisation. 

To first order the reionisation process can be modelled by 
equating the number of photons produced per baryon to the number 
of recombinations in the ionised IGM (e.g. Madau et al. 1999J. The 
volume fraction of ionised hydrogen, Qn u is then given by 



dQmi 
dt 



N ion 

«H,0 



?HiiCn H ,oa(T)(l + ^) 3 , 
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where JVi on is the total number of ionising photons available for 
reionisation per second per comoving Mpc, nn.o = 1.90641 x 
10~ 7 cm~ 3 is the current mean number density of hydrogen in the 
IGM, a(T) is the recombination coefficient of hydrogen, which is 
a function of the IGM temperature T, C = {fig)/{nH) 2 is the 
clumping factor of the gas in the IGM and z is the redshift. We as- 
sume that the ionised gas in the IGM has a temperature of 20,000 
K, while we compute the clumping factor of the IGM gas from the 
simulation, using only gas with overdensity A < 100, thus ex- 
cluding gas inside dark matter haloes (because recombinations in 
that gas are already taken into account in the escape fraction cal- 
culations). We find a redshift-dependent clumping factor between 
1.5 and 6.5. consistent with previous studies (Pawlik et al. 2009; 
|Shull et al.|2 012 ). We compute the Thomson optical depth, which is 
the quantity measured by the WMAP satellite, by integrating Qh ii 
over all redshifts: 
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with <jt the cross section for Thomson scattering and c the speed 
of light. 



3 THE CONTRIBUTION OF PROTO-GALAXIES TO 
COSMIC REIONISATION 

In Fig. [T] we show the average escape fraction as function of red- 
shift for proto-galaxies in haloes of different masses. The average 
escape fraction rises with time, but proto-galaxies inside a certain 
mass halo at the same redshift may have very different escape frac- 
tions, for example due to a different formation history or environ- 
ment, which causes the large standard deviation in the mean. Due to 
the efficiency of stellar feedback in clearing away the gas from the 
dense sites of star formation, the escape fraction in 10 7 Mq haloes 
is higher than in the 10 8 Mq haloes. Proto-galaxies in haloes with 
masses above 10 9 Mq have, at all redshifts, significantly lower es- 
cape fractions than their counterparts in lower mass haloes, due to 
their larger and denser gas content. 

Haloes with masses below 10 s Mq dominate the ionising pho- 
ton budget at redshifts higher than 10, due to their high escape frac- 
tions and high abundance. After redshift 10, the ionising emissivity 
of the proto-galaxies in these haloes drops as a result of suppression 
of star formation by the uniform UV-background. The background 
heating only suppresses star formation in haloes that do not contain 
enough dense gas to shield against the radiation. This counteracts 
the effect of the high escape fraction, since ionising radiation is 
mainly produced by massive, young stars and suppression of star 
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Figure 1. The escape fraction of ionising photons as function of redshift. 
The black solid line represents the escape fraction averaged over all haloes, 
while the blue dashed, green dotted and red dash-dot lines denote the es- 
cape fraction averaged over proto-galaxies in haloes with virial masses of 
10 7 Mq, 10 8 Mq and 10 9 Mq, respectively. The grey area represents the 
standard deviation of the mean. 




z 

Figure 2. The contribution of proto-galaxies to reionisation. In all panels 
the black solid lines represent all haloes, while the blue dashed and green 
dotted lines show the contribution of haloes with masses below 10 7 M 
and 10 8 Mq, respectively. The red dash-dot lines represent the contribution 
of Pop III stars. Top panel: The halo mass below which 50% of the ionising 
photons is produced, Mo. 5, at each redshift. The grey area represents the 
range of halo masses in which star formation is taking place. Middle panel: 
The optical depth for Thomson scattering, r c as function of redshift. The 
data point represents the value of T e as found by measurements with the 
WMAP satellite I Komatsu et al. 201 1 1. Bottom panel: The volume filling 
fraction of ionised hydrogen, Qh II , as function of redshift. 

formation results in little or no ionising radiation being produced 
(see Fig. [2j. Low-mass haloes therefore only contribute to the ion- 
ising photon budget when the part of the Universe they live in is 
still neutral, while more massive haloes have enough dense gas to 
shield against the external radiation and continue to form stars. 

In Fig. [2] we show the evolution of Qnu and r c with redshift. 
Reionisation is complete at redshift 10.5, with a duration (defined 
as the redshift interval in which Qu 11 changes from 20% to 80 %) 
of Az ~ 2.1. This is before the lower limit of the end of reionisa- 




Figure 3. The cumulative number of ionising photons per baryon as func- 
tion of redshift. Different colors denote the contribution of proto-galaxies 
below a certain stellar mass. In our computation full reionisation requires 
1.6 photons per baryon, which is represented by the black dotted line. 



tion set by the measurements of quasar spectra. The dashed curve 
shows the reionisation history including only proto-galaxies resid- 
ing in < 10 7 Mq haloes, which do not produce enough photons to 
reionise the Universe. The bulk of photons is produced by proto- 
galaxies in haloes with masses between 10 7 Mq and 1O 8 M0. 
Reionisation is only delayed by Az ~ 0.5 if we exclude all sources 
in haloes with masses larger than 10 8 Mq. 

The integrated Thomson scattering optical depth that we find 
in our model is r c = 0.096, which is well within the error bars 
from the WMAP measurement (Kom atsu et al.|2011| l. This shows 
that the reionisation history we find is consistent with the two main 
observational constraints of reionisation, the absorption features in 
high-redshift quasars and the Thomson optical depth as observed 
by WMAP. 

The contribution of Pop III stars to reionisation is negligible in 
our simulations. Although these sources produce copious amounts 
of ionising photons (Schaerer 2002), they are short-lived and not 
abundant enough to contribute significantly to reionisation. The 
contribution of Pop III stars to the total photon budget is exceeded 
by metal-enriched Pop II stars at redshifts below 15. We thus con- 
clude that although reionisation started with the appearance of the 
first Pop III stars, they did not contribute significantly to reionisa- 
tion on global scale. 

In the top panel in Fig. [2] we show the halo mass below which 
50% of the ionising photons are produced as function of redshift. 
At redshifts higher than 10, half of the ionising photons are pro- 
duced by proto-galaxies in haloes with masses between 10 7 Mq 
and 10 s Mq. After this redshift, haloes with higher masses take 
over the photon production, because star formation is suppressed 
in the lower mass haloes due to the UV-background. 

To get a better picture of the mass of galaxies mostly con- 
tributing to reionisation, we show in Fig. [3] the cumulative num- 
ber of photons per baryon produced by proto-galaxies below a cer- 
tain stellar mass as function of redshift. Given our estimates of the 
ionising emissivity and the clumping factor, reionisation requires 
approximately 1.6 photons per baryon. Proto-galaxies with stellar 
masses below 10 6 Mq have produced this number of photons by 
redshift 10, with most ionising photons being produced by proto- 
salaxies with stellar masses between 10 5 M Q and 10 6 M . 
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Figure 4. The observability of the sources of reionisation in the UV contin- 
uum and three recombination lines. In all plots colors are similar to the col- 
ors in Fig. [3] The shaded area denotes the maximum and minimum bright- 
ness and the orange area covers the range detectible by the relevant instru- 
ment on JWST. 



4 THE OBSERVABILITY OF THE SOURCES OF 
REIONISATION 

One of the science goals of the James Webb Space Telescope 
(JWST) is to study the galaxy population during reionisation JGard-| 
|ner et al.|2006| l. The deepest NIRCam survey will search for high- 
redshift objects with exposure times of 10 6 seconds, resulting in a 
flux limit of 1.4 njy at 2 (jm. This can be converted into AB mag- 
nitude using 



mm = 31.4- 2.5 log(/„ 



(3) 



resulting in a magnitude limit of ttiab = 31.0. In addition, JWST 
will search for recombination line radiation from the first galaxies 
using the MIRI and NIRspec instruments. At 5.6 /im, the limiting 
flux of the MIRI instrument is 23nJy or diab = 28, while NIRSpec 
will measure line intensities down to 2 x 10 -19 erg cm -2 s _1 with 
a resolution 7? = 1000. 

We determine the flux that the NIRCam survey would receive 
from the proto-galaxies in our simulation by computing the spec- 
tra of the Pop II and Pop III stars for radiation with energy be- 
low the ionisation energy of hydrogen. In Fig. [4] we show the red- 
shifted flux from the sources in the 2 fim band, although due to 
the flatness of the spectrum in this wavelength range the results do 
not depend sensitively on the choice of wavelength. Reionisation 
is driven by proto-galaxies with stellar masses around 1O 6 M0. 
At redshift 10.5, when the proto-galaxies have produced enough 
photons to reionise the volume, these sources are not brighter 
than f v = 0.09 njy or jtiab = 34 , making it impossible for 
JWST to observe them. If we rescale the flux limits according to 
/ii m oc 1/y/texp, exposure times of 2 x 10 8 s are necessary to ob- 
serve these sources. Since we do not account for dust attenuation 
between the proto-galaxy and the observer, we can only give upper 
limits of the flux. 

In addition, we show in Fig. [4] the brightness of the proto- 
galaxies in our simulation in three emission lines, Ha, Ly-a and 
Hell 1640A. The flux in these lines is given by 1 Johnson et al. 
|2009l > 

X crn (l + z)R 



/(Aob 



4ncDl(z) 



(4) 



where £ eia is the luminosity along the line of sight, Df j (z) is the lu- 
minosity distance at redshift z, c is the light speed and R = A/ A A 
is the spectral resolution. Assuming that the galaxy is unresolved, 
we compute l cnl for every proto-galaxy by summing up the contri- 
bution to the emissivity from every gas particle. 

JWST will observe the Ha line with the MIRI instrument. 
At redshift 12 and below, proto-galaxies of M* ^ 10 6 M© are 
bright enough to be detected with MIRI. The Ly-a and Hell 1640A 
lines will be observed with the NIRSpec instrument. The greater 
sensitivity of the NIRSpec instrument could in principle make it 
possible to observe the Ly-a emission line from all sources with 
Mi, ^ 10 4 Mq at redshift 20 and below. However, this line is 
difficult to detect at these high redshifts due to the large Gunn- 
Peterson optical depth. Scattering off interstellar neutral gas could 
make the line observ able even through a fully neutral IGM (Dijk- 
stra & Wyithe||2010| >. The flux in the Hell 1640A line is always 
lower than the Ha flux. It is therefore most likely that the sources 
of reionisation will be observed in the Ha or Ly-a line. 



5 CONCLUSIONS AND DISCUSSION 

We have presented high-resolution cosmological simulations of 
galaxy formation, from which we calculated the ionising photon 
production and escape fraction of a large statistical sample of galax- 
ies. From these simulations we computed the contribution of proto- 
galaxies to cosmic reionisation in a self-consistent way. Our main 
findings are: 

• Reionisation is primarily driven by proto-galaxies in dark- 
matter haloes with masses between 1O 7 M0 and 10 s Mq, which 
have very high escape fractions, because supernova feedback effi- 
ciently acts to clear away gas from the sites of star formation. 

• Star formation in these haloes is suppressed by UV-feedback, 
these proto-galaxies therefore only contribute to reionisation when 
the part of the Universe they live in is still neutral. 

• After reionisation the Universe is kept ionised by massive 
galaxies with lower escape fractions, that start appearing more fre- 
quently. 

• We find a strong mass and redshift dependence of the escape 
fraction that suggests that galaxies above the observational limits 
of present surveys do not contribute enough photons to drive reion- 
isation. 

• Pop III stars do not contribute significantly to the reionisation 
process. 

• There is great prospect that JWST will observe the proto- 
galaxies that reionised the Universe in the Ha and Ly-a recom- 
bination lines. 

In our galaxy formation simulation radiative feedback on the halo 
gas from sources within the halo is neglected. Radiative feedback 
is capable of evacuating gas from haloes with Af v ; r < 10 7 Mq 
( |Wise & Cen|2009| >, thereby suppressing further star formation in 
the halo. We find that supernova feedback has the same effect, albeit 
with a short delay. Since the time it takes for the dense, cold gas in 
the halo to be converted into stars (~ 1 Gyr) is much larger than 
the lifetime of the massive stars that end their lives as supernova 
(~ 1 Myr for pair-instability supernovae), this delay will not affect 
our results significantly. 

We model reionisation with a uniform UV-background, dis- 
regarding the contribution from local sources. We could therefore 
be underestimating the suppression of star formation in haloes that 
are not able to shield against the ionising radiation. This likely 
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occurs in haloes with masses Mvii < 10 Mq i Okamoto et al. 



[2008) Hasegawa & Semelin 2012). Since this will only affect 
proto-galaxies in a region of the Universe that has already been 
ionised, this will not change our conclusions. However, it could 
delay the completion of reionisation due to suppression of clus- 
tered low-mass sources and reduce the number of photons avail- 
able for reionisation. Assuming all star formation in haloes with 
-Afvir < 10 9 Mq is suppressed after reionisation brings the cumu- 
lative number of ionising photons per baryon at redshift 6 in our 
simulation close to the observed value (Bolton & Haehnelt 2007). 

The resolution in our simulation is around 6 physical pc at red- 
shift 15, which is high enough to resolve regular and giant molecu- 
lar clouds at the sites of star formation. However, we do not resolve 
the birth cloud of the stellar population, which would lower the es- 
cape fraction. Since this equally applies to all haloes in the volume, 
it does not affect our conclusion on the mass range that contributes 
most photons to reionisation. The effect of a lower escape fraction 
can be estimated by scaling the curves in Fig. [^accordingly. 

The volume of our simulation is 4 Mpc, which could possibly 
bias our results. The initial conditions for the cosmological simula- 
tion were chosen to avoid biased regions with high a peaks. This is 
reflected in the mass function of dark matter haloes in our volume, 
which is in agreement with the average of the Sheth-Tormen mass 
function ( |Sheth & Torm en 2002). Although we cannot accurately 
sample the high-mass end of the halo mass function, the strong 
mass-dependence we find for the escape fraction indicates that the 
contribution from these rare high-mass sources is small. 

Our results show that reionisation is initially a local process 
driven by the appearance of many low-mass haloes which do not 
cluster strongly and ionise their neighbourhood until H II regions 
overlap and cover larger cosmic volumes. 
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